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Abstract 

A multi-channel generalization of Doniach's Kondo necklace model 
is formulated, and its phase diagram studied in the mean-field approxi- 
mation. Our intention is to introduce the possible simplest model which 
displays some of the features expected from the overscreened Kondo lat- 
tice. The N conduction electron channels are represented by N sets of 
pseudospins Tj, j = 1,...,N, which are all antiferromagnetically coupled 
to a periodic array of \S\ = 1/2 spins. Exploiting permutation symme- 
try in the channel index j allows us to write down the self-consistency 
equation for general N. For N > 2, we find that the critical tempera- 
ture is rising with increasing Kondo interaction; we interpret this effect 
by pointing out that the Kondo coupling creates the composite pseu- 
dospin objects which undergo an ordering transition. The relevance 
of our findings to the underlying fermionic multi-channel problem is 
discussed. 
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1 INTRODUCTION 



The essential physics of mixed valent and heavy fermionic systems jl| is habit- 
ually described by some suitable version of the periodic Anderson model ||. 
In this framework, heavy fermionic systems correspond to the limiting case 
of nearly integral valence, which permits us to use a Kondo lattice model for 
the approximate handling of the low-energy, low-temperature behaviour of 
the Anderson lattice [pi. Though already several steps removed from physical 
reality, it is widely accepted that the Kondo lattice model incorporates the ba- 
sic ingredients necessary to understand the competition between two opposing 
tendencies: the formation and eventual ordering of magnetic moments, and the 
formation of a heavy Fermi sea. Research has for a long time been focussed on 
the case where localized spins S — 1/2 are compensated by electrons moving 
in a non-degenerate band 

H KL = J2 <k)c+ CT c kff + ± £ £ e^~^ £ c+^c^S (1) 

k,cr Zlj g k,k' a/3 

where the first term describes a non-degenerate band, and the second the 
antiferromagnetic (J > 0) Kondo coupling at each site g. L is the number 
of lattice sites, and cr is the vector of Pauli matrices. H^l is the simplest 
version of the Kondo lattice hamiltonian; we expect that it shows the same 
overall features as the S = N/2 models, where N is the number of conduction 
electron channels participating in the Kondo screening process. This special 
case of exact compensation seems to be the most favourable for the formation 
of an overall singlet, i.e., a heavy Fermi liquid without any residue of magnetic 
order. Even these perfectly screened models are, of course, capable of magnetic 
ordering 0, [|, |, || , but the arising of a completely non-magnetic state can at 
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least be claimed to be a natural option 0. 

In recent years, it has become increasingly clear that considering the un- 
derscreened S > N/2, and overscreened S < N/2, models brings not just 
additional complications but interesting new physics, and it is also essential to 
describe many, if not all, of the actual heavy fermionic materials. 

The properties of an isolated underscreened, or overscreened, Kondo im- 
purity || have been studied in great detail, and for a number of physical 
quantities, exact results are known ||, [UJ. Still, understanding can not yet be 
considered complete; in particular, we could wish for a clearer physical picture 
of the low-temperature behaviour of the overscreened model with its intriguing 



non-integer ground state degeneracy JTT], |12 |. 

Much less is known about magnetic ordering in the corresponding periodic 
models. Just as in the case of perfect screening S = N/2, there must be 
a competition between Kondo screening and intersite interactions; however, 
the picture is complicated by the fact that the local object which tends to 
arise from the Kondo effect, has a more intricate internal structure. For the 
underscreened case, even a very strong Kondo coupling can not do better than 
to leave a reduced spin S — N/2; thus one is led to consider the ordering of the 
residual magnetic moments. The underscreened Kondo lattice is apparently 
destined to become magnetic and it was suggested that its study should 
lead to understanding heavy fermion magnetism in general [ITJ . 

The situation is much less clear in the overscreened case. The model which 
was first introduced as a matter of academic interest ||, was later argued to 
represent the physical situation for two-level systems interacting with several 
conduction electron channels [15, HJ, and for U-based heavy fermion systems 
||17|| . Subsequent experiments on structurally disordered metallic nanoconstric- 
tions and on the Kondo alloy Yi_ x XJ x Pd 3 |TP|], gave ample confirmation 
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of the applicability of the overscreened Kondo model. Another remarkable 
realization of the multi-channel Kondo effect is provided by the transport 



properties of Pbi^Ge^Te ||2(| . — The solution of the single-impurity problem 
shows a puzzling mixture of features which one would naively associate either 
with the absence, or the completion, of spin compensation. The ground state 
possesses a zero point entropy [|1(], [11], |T^| which, in the limit of an infinite 
Fermi sea, does not correspond to a half-integer spin, but reduces to In (2) if 
the system is finite. The finding of a zero point entropy might lead us to think 
that the ground state is possessing a residual magnetic moment; however, cal- 
culations show [§, |l(| that magnetism is quenched as T— >0, albeit much more 
slowly than in the S = N/2 case. All in all, it turns out to be fairly difficult to 
form an intuitive picture of the multichannel Kondo system whose behaviour is 
governed by an intermediate-coupling fixed point. Correspondingly, we have 
only the haziest idea of what a periodic array of overscreened Kondo centres 
is supposed to be doing. 

In either case, when we undertake the study of the magnetic phases of the 
periodic models, we have to deal with the ordering of composite spins created 
by a local Kondo coupling; and the degree to which Kondo compensation can 
progress is itself limited by the intervention of intersite interactions. Curiously, 
the effect can go opposite ways in the underscreened, and overscreened, cases. 
For the underscreened lattice, the resulting total spin is necessarily less than S, 
so the ordering temperature is found to show a decreasing tendency with 
increasing J. In the overscreened case, the Kondo effect is actually building 
up a screening object which can be ultimately larger than the screened spin 
S. We are going to find this tendency for sufficiently large (N > 3) numbers 
of channels. 
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2 PSEUDOSPIN MODELS: THE KONDO NECKLACE AND ITS GEN- 
ERALIZATIONS 



The fact that the conduction electrons are forming a Fermi sea, with a large 
number of arbitrarily low-lying excitations, is essential for the appearance of a 
non-analytic energy scale [the Kondo temperature T K oc Bexp (— 1/Jp(ep)), 
where B is of the order of the bandwidth, and p{ef) the density of states at 
the Fermi level] in the theory of the single-ion Kondo effect. It is, as yet, an 
open question whether there is a non-analytic energy scale associated with 
the formation of a non-magnetic ground state of the periodic Kondo model 
(or, to put it more pointedly, whether there is a Kondo-effect in the Kondo 
lattice and if yes, whether the "lattice Kondo temperature" is different 



from the single-ion Tk f| 0, ^31 - It can be, however, argued [24] that the 



competition between the Kondo effect, and the interactions opposing it, can 
be successfully mimicked by just considering the gross effects of different cou- 
plings: one can have a crude version of the phase diagram of the Kondo lattice 
without the "true" Kondo effect. One can then forget about the low-lying 
electron-hole excitations which can lead to infrared divergencies, and consider 
the conduction electrons only to the extent that they provide spins which tend 
to be aligned antiparallel to the localized spins. One can thus reduce the mixed 
spin-fermion problem to a pure spin problem. This was the rationale behind 



introducing Doniach's [25] necklace model 



H KN = JJ2 S m ■ r m + W £ + 7* 7*) (2) 

m (mn) 

where the pseudospins |r| = 1/2 represent the spin degrees of freedom of the 
conduction electrons, and W is a characteristic amplitude of their propagation. 
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Hkn is the pseudospin version of the S = 1/2 lattice model given in eqn. (1); 
clearly, W is an effective parameter which depends in a complicated manner 
on the parameters of the underlying electronic model. The mean field treat- 
ment of Hkn by Doniach J2f| led to a ground state phase diagram which looks 
qualitatively the same as expected from a (still missing) complete treatment of 
(1), for the case of a half-filled conduction band: the ground state is antiferro- 
magnetic for J/W < 1, and non-magnetic (apparently Kondo-compensated) 
for J/W > 1. We conclude that the study of necklace-type pseudospin mod- 
els is a useful prelude to the investigation of the full fermionic Kondo lattice 



problems. — We note that recently, Strong and Millis p6| investigated a gen- 
eralized form of the 5 = 1/2, N — 1 necklace model with the ultimate purpose 
of gaining insight into the behaviour of heavy fermionic systems. 

We should emphasize that we use the term "necklace" , as opposed to "lat- 
tice", models to signify that the conduction electron sea has been replaced 
by a set of pseudospins, and not to mean that the model is necessarily one- 
dimensional. In (2), the sum over (mn) is over nearest-neighbour (nn) pairs 
in any lattice. In fact, usually we will resort to the mean field approxima- 
tion (MFA), so lattice dimensionality plays no essential role. The results 
obtained in the MFA are thought to be a reasonable approximation for the 
three-dimensional models, but they should serve as a rough guide even in one- 
dimension P7| , except that whenever long-range order is found, it should be 
interpreted as indication of quasi-long-range order. 

Considering the rudimentary stage of the understanding of the under- 
screened, and quite particularly, the overscreened, periodic Kondo models, we 
can appreciate the need for any insight to be gained from the study of related 
simpler models. In ||28| , we introduced the under screened, spin-S 1 necklace 
models, and discussed their properties. The underscreened necklace models 
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turn out to be inherently antiferromagnetic, with an easy-plane anisotropy; 
the Kondo coupling merely influences the size of the ordered moment. The 
relative ease of the solution is thanks to the fact that the dimensionality of the 
Hilbert space is merely oc S. 

In the present work, we perform a similar investigation of a class of over- 
screened necklace models. We introduce the hamiltonian 

j N w N 

H = k £ 5m ■ £ ^ + T7 E £(^ m T*„ + rj>jg (3) 

m 3=1 (mn)j=l 

where iV different kinds of pseudospins Tj for j = 1, ...,N are meant to repre- 
sent the N screening channels. The localized moments (in truth, /-spins) are 
taken to be \S\ = 1/2. The x — y coupling which should correspond to the 
propagating character of the conduction electron spin degrees of freedom, acts 
independently in each channel. In order to have a meaningful large-iV limit 
of the model, one has to keep the individual couplings of order 1/N. We are 
restricting ourselves to the case when all channels are equivalent; one should 
be aware, though, that models with inequivalent channels |30| should be of 



great physical interest. The study of some such models is in progress. 

A general comment about the necklace replacements of the true Kondo 
problems should be made here: in the necklace models, the Kondo effect, and 
the Kondo bound state, become very local. This corresponds to the actual 
physical situation for very strong Kondo couplings. Thus, the necklace models 
can be used in good faith for S>k/2, where the large-J limit is continuously 
connected to the regime of small J's. The finding of an intermediate coupling 
fixed point for S < k/2 tells us that this is assuredly not the case for an 
overscreened Kondo impurity. However, we would still like to argue that the 
study of the lattice model (3) should have relevance for the overscreened peri- 
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odic Kondo model. It can be generally argued that in Kondo lattice models, 
especially in the presence of ordering, intersite effects may prevent us from 
reaching the fixed point of the corresponding single-site models, and so the 
subtleties of the impurity solution may become irrelevant for the lattice. This 
argument should apply quite particularly to the overscreened case where the 
behaviour is known to be extremely sensitive to the presence of external fields 



T0|, |29|. It has been pointed out by Andraka and Tsvelik pTJ that intersite 



interactions, acting like fields, can be expected to turn the system towards 
a strong-coupling fixed point. An important piece of corroborative evidence 
is that for the two-impurity problem, destabilization of the "marginal Fermi 
liquid" single-impurity behaviour was found by Ingersent, Jones and Wilkins 
f32fl . Thus we think it not implausible that in periodic models, ordering of 
tightly bound objects created by the multichannel Kondo coupling can be tak- 
ing place, and then our model hamiltonian (3) acquires relevance. 



3 MEAN-FIELD THEORY OF THE OVERSCREENED NECKLACE 
MODEL 

The overscreened model presents much larger difficulties than the underscreened 
one 



28| because the dimensionality of the Hilbert space is now exponentially 
large in N. We are, however, aided by the huge symmetry of H: neglecting 
lattice translations, the symmetry group is G = SO(3)<2)Sn- H has full spin- 
rotational symmetry, as well as invariance under arbitrary permutations of the 
channel indices j. We will find that symmetry simplifies the appearance of the 
spectrum to such an extent that at least the single-site mean field problem 
becomes solvable for general N, without any further approximations. 

Assuming a simple magnetic order in which all neighbours of a given site 



S 



have their r-spins polarized in the same direction, we can do the familiar MF 
decoupling to arrive at the single-site hamiltonian 




N 



N 



N 



(4) 



where u = Wz(t) is the mean field strength, z being the coordination number. 
For the sake of solving just the MF problem, we rotated the quantization axis 
to align up with the mean field; in the sense of the full hamiltonian H, this 
must be an arbitrary direction in the x — y plane. For bipartite lattices, there 
is no formal difference between ferromagnetic W < 0, and antiferromagnetic 
W > 0, intersite couplings; to have a closer correspondence with physical 
reality, we will have to take W > 0. 

The Kondo term in H MF has the single-site version of the high symme- 
try described before: it is invariant under arbitrary rotations of the total spin 
r] = S + J2f=i r j) as we ll under permutations of the N channels. Switching 
on the effective field destroys spin-rotational invariance, but we will still find 
it useful to think of the states as derived from splitting the highly degenerate 
levels of the underlying S'0(3)-invariant problem. In any case, we have still 
got the permutational symmetry. Actually, this rests on having chosen a j- 
independent effective field: 



In principle, MF solutions breaking the Sn symmetry are imaginable. For small 
values of N, we were exploring this possibility, and found that the symmetrical 
solution is more advantageous, so we feel confident in using the form of Hmf 




(r z N ) = (r) 



(5) 
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specified in (4). 

In what follows, we set out to classify the eigenstates according to the 



irreducible representations of Sn, by constructing their Young-tableaux [33 



As it turns out, this classification is intimately related to the classification 
according to the total spin r], and its ^-component r/ z , so it is useful to start 
with the high-symmetry problem u = 0. 

Let us start with the unique rf = (N + l)/2 state 

N 



(6) 



where ft denotes the 5-spin, and the r-spins enter the Young tableau (YT) for 
the identity representation of Sn- Obviously, the S'-spin cannot be included 
in the YT because there is no symmetry with respect to interchanging S with 
one of the rs. Following custom, we will refer to the r-part of such a graphical 
representation as a YT, while the whole picture will be called a "diagram"; 
since with the inclusion of the S'-spin, it is not a YT in the conventional sense. 

(6) is the r] z = (N + l)/2 component of the maximum spin rj = (N + l)/2 
multiplet. Its energy is 



In the subspace = (N — l)/2 we are going to meet several kinds of 
states. The permutationally symmetrical ones are 
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N 



(8) 



and 



N-l 



it 



( 






T 




T 


1 



(9) 



As usual, YT represent states fully symmetrized along a row, so it is unambigu- 
ous what (8), and (9), stand for. These are not eigenstates of the total Kondo 
coupling; the actual eigenstates are their properly chosen linear combinations. 
One of these must be the if = (N — l)/2 component of the rj = (N + l)/2 
multiplet. The other belongs to a symmetrical rf = (N — l)/2 multiplet, with 
the energy 

N + 2 

eo = (10) 



This is actually the ground state energy of the local Kondo coupling. The 
corresponding eigenstate can be roughly described by saying that the S-spin 
is pointing antiparallel to an assembly of r-spins. Actually, it is not quite 
that: it contains an admixture of states (9) in which a r-spin is antiparallel 
to all other spins. 

The ground state of the Kondo-term is iV-fold degenerate. We can inter- 
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pret this as a result of an effective ferromagnetic coupling which the common 
Kondo-coupling to the S-spin induces between the r-spins. This feature is 
not unlike to that found in rigorous treatments of the overscreened Kondo im- 
purity jlOj ^9[. — Note, however, that here it leads to ascribing the zero-point 
entropy /c^lniV to a Kondo site and thus , as discussed in the Introduction, 
we fail to recover the subtle features of the single-ion solution of the origi- 
nal overscreened Kondo problem. However, we have the correct features for a 
strong-coupling solution which may turn out to be of relevance for the lattice 
case. 

€q and 6q can be seen to result from coupling the S = 1/2 spin either 
parallel, or antiparallel, with a fully polarized T = |X)j T j| = N/2 set of 
r-spins. 

All the other N — 1 states in the rf = (N — l)/2 subspace belong to the 
only other Young tableau which can be constructed with one of the r-spins 
down: 







N-l 














T 




T 




I 





(11) 



Antisymmetrizing along the column is effectively the same as putting two r- 
spins in a singlet combination, permitting to line up only the residual N — 1 
spins. The dimensionality of the representation specified by the YT above is 
just d\ = N — 1, accounting for all the remaining = (N — l)/2 states which 
are thus found to be degenerate. Here we first see the enormous simplification 
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brought by the S^-symmetry: spin-rotational symmetry alone would permit 
the existence of iV different levels corresponding to rj = (N — l)/2 multiplets; 
SAT-invariance tells us that there are just two levels, corresponding to the two 
different possible shapes of Young tableaux. One of these was inherited from 
the rf = (N + l)/2 subspace and corresponds to T = TV/2; the other is the 
one in (11), with a down-spin in the second row, for which T = (TV — 2)/2. 

- For our particular problem, which has to do with permutational symmetry 
acting within 577(2), all YT have at most two rows )[33|| . 

It is crucial to observe that this scheme holds generally: each time when 
we step down r] z , only one new YT appears, namely the one where one box is 
removed from the upper row, and a down-spin box is added to the second row. 
This corresponds to a new value of the total r-spin T, which is 1 less than 
the previous lowest T-value. Though T = J2j T j is n °t a conserved quantity 
for the Kondo-term, definite values of T are unambiguously identified with 
definite shapes of Young tableaux, T = (TV — 2m)/ 2 corresponding to a YT 
with N — m boxes in the upper row. 

For vf = (TV — 2m + 1)/2, most of of the states can be represented by turn- 
ing down a spin in the first row of a YT which we had seen previously: these 
are stepped-down versions of rj > if states, corresponding to levels which had 
been identified previously. The only new diagrams are 



N-m+1 



T 




T 


I 




I 





m-l 



(12) 
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associated with combining T = (N — 2m + 2)/2 antiparallel with S — 1/2, i.e., 
with r] — (N — 2m + l)/2, and the diagram with the new YT 

N-m 

i 1 



ft 



T 




T 


I 




I 





V m ' (13) 

describing the parallel alignment of T = (N — 2m)/2 total r-spin with the 
S-spin. This also has 77 = (iV — 2m + 1)/2, but belongs to a new representation 
of Sjvj along with the new value of T. The diagram is uniquely associated with 
the maximum-?? 2 component of the newly found multiplets which are therefore 
all degenerate. The degeneracy of the T = (N — 2m)/2, rf = (N — 2m + l)/2 
level is 

dm = Cm - Cm-1 = Cm +1 JV~+1 

where the Cs denote binomials. 

The energies are simply given by 

JT-S=±[7i( V + l)-l-T(T + l)] (15) 
The degeneracy of each energy level is the product of the corresponding d m , 
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and the S'0(3)-associated multiplicity 2rj + 1. 

For N even, the procedure finishes with 

N/2 

i ^ 



T 




T 


I 




I 



(16) 

and a similar diagram with the .S-spin JJ*. These stand for a rf^-dimensional 
subspace of r-singlets, leaving the S'-spin free to make a disconnected doublet, 
with energy e^ /2 = 0. — For N odd, the last diagram is 



(N+l)/2 



I 






T 




T 


T 


1 




1 






i 






(N-l)/2 







signifying an overall singlet r] = 0. 

In either case, two energy levels belong to each form of the YT with < 
m < [(N + l)/2] - 1, and one level goes with m — [(N + l)/2]. The total 
number of energy levels is N + 1. 

In the presence of an effective field w^O, the symmetry is reduced to Sn, and 
i] is no longer a good quantum number. However, the full spectrum can still 
be determined, because the eigenvalue problem is easily seen to separate into 
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a number of of two-, and one-dimensional problems, essentially because the 
S = 1/2-spin can be stepped only once. Permutation symmetry requires that 
only states with the same shape of the YT mix. Easiest is the case m — 0, 
with a YT consisting of a single row of | _ , and j -spins. In the orthonormal 
basis 

N-r r 



r 



T 




T 


1 




1 



E 



jl<...<jr 



T ••• I ••• I ••• T) 



(18) 



and 



N-r+1 



r-1 



T 




T 


I 




I 



y 

/r<N . ^ . 

V U r-1 Ji<-<>- 1 



T ••• I ••• I ••• T) 



-£?mf is represented by the matrix 

( (1/AN)(J -2to)(N -2r) ( J/2N) y/r(N - r + 1) 



(19) 



(20) 



V (J/2N)Jr(N -r + 1) -(1/4N)( J + 2u>)[N - 2(r - 1)] J 



The ground state is now non-degenerate, corresponding to r = 1, and a max- 
imum polarization along the effective field. Note that while, in the absence of 
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an external field term, the pseudospin model gives only a very poor imititation 
of the highly non-trivial quantum-mechanical ground state of the electronic 
overscreened model, we have much less reason to doubt the validity of the pic- 
ture obtained from the necklace model for moderately strong external fields. 
Since in a lattice, intersite interactions amount to a field acting at any par- 
ticular Kondo site, herein lies our hope that the multichannel necklace model 
can give some guidance as to the behaviour of the multichannel Kondo lattice. 



We have already quoted arguments pT|, showing that intersite interactions 



can make the system turn towards a strong-coupling fixed point: this kind of 
behaviour can be imitated by a necklace-type model. 

Similarly, for states with m .[-spins in the second row of the YT, the two- 
dimensional subspace is spanned by 



N-m-r 



T 




T 


1 




I 


I 




I 





/// 



(21) 



and 
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N-m-r+1 



r-1 



T 




T 


I 




1 


1 




1 





m 



(22) 



The eigenvalue problem is similar to that of m — 0, only we have to remem- 
ber that in a YT, the states are antisymmetrized along the columns, so the 
symmetrization of |- and |-spins along the first row is not allowed to bring 
J, -spins into the first m boxes. Effectively, in (20), in the off-diagonal elements 
N must be replaced by iV — 2m, giving the matrix 



/ (l/4N)(J -2lu)[N -2(m + r)} (J/2N)^Jr(N - 2m - r + 1) \ 

v (J/2N)^r(N -2m-r + 1) -(l/4iV)(J + 2u)[N - 2(m + r - 1)] j 



(23) 



The solution of the eigenvalue problem would be trivial to write down, but we 
do not need the lengthy expressions. First of all, we are interested in the mean 
field result for the magnetic transition temperature T N , for which linearized 
eigenvalues 



A+ - J (N 2m) (N-2m)(N + l-2m-2r) 
A ™ " 4N {N ~ 2m) 2N(N + 1 - 2m) U 



(24) 



and 



18 



J , n n , (iV + 2-2m)(iV + l-2m-2r) 
A"„ = - — (iV + 2 - 2m) - ^ — ^ ^ (25) 



AN 



2N(N + 1 - 2m) 



are sufficient. We also need (r) for the corresponding eigenstates 



+ _ (N-2m)(N + 1 -2m-2r) 4rw iV + 1 - 2m - r 
Tm*. — ttz : ~ ; t~t~: : ; rr~ l^'-v 



rnr 



2N(N + 1 - 2m) 



J7V (iV + l-2m) E 



and 



(iV + 2-2m)(iV + l-2m-2r) 4ru N + 1 - 2m - r 

T = h 



2iV(iV + 1 - 2m) 



JiV (iV + l-2m) 3 



(27) 



In the linearized self-consistency equation for (r), the denominator is just the 
uj = value of the single-site partition function 



Z = Yl d m (N - 2m + 2)e- ,{N - 2m V JNT " + (N - 2m )e^ N - 2m+2 ^ ANT - 



m=0 



(28) 



while the numerator is the suitably weighted sum of (26) and (27) 



1 2 

0") = ^" S rfm 
Z m=0 



N-2m N-2m+l 
, r=l r=0 



(29) 



Actually, a slight complication has to be dealt with before arriving at (29). In 
the second r-sum, the terms 1 < r < N — 2m arise from the two-dimensional 
eigenvalue problem (23). However, the state with r = 0, and S-spin ff, is not 
coupled to any other state, it presents a one-dimensional eigenvalue problem. 
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Similarly for r = N — 2m + 1, and the S'-spin JJ-. These are actually the 
f] z = ±(iV — 2m + l)/2 components of the corresponding rj = (N — 2m + l)/2 
multiplets, and are thus eigenstates by construction. It is, however, easy to 
check that these cases are accounted for by extending the second sum in (29) 
to include r = 0, and r = N — 2m + 1. 

Though the r-sums in (29) are easily done (with the result that (r) disap- 
pears from the equation), the m-summation is not, and we did not find any 
transparent from into which (29) could be cast. Essential simplification can be 
achieved only for the limiting case J/W — > oo when it is sufficient to consider 
the action of Hmf within the subspace of the rj = (N — l)/2 Kondo ground 
states associated with the diagrams shown in (8), and (9). From (29), we find 

Tn= W(N- l)(N + 2? 
j/w^oo N 6 N 2 (N + 1) V ; 



The above result can be easily interpreted by noticing that within the 
Kondo ground state set of states Tj acts like 

N + 2 

T > ivWTT) 77 (31) 

independently of j. Thus, in the restricted Hilbert space formed by taking 
the direct product of the Kondo ground state sets for each site, the intersite 
coupling term in (3) can be exactly represented by the effective hamiltonian 

heff = W ^ N + ^ )2 E m + « (32) 
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(30) is just the mean-field solution for (32). Note that it gives for N 



this is in accordance with Doniach's finding |25| of a transition to the fully 
Kondo-compensated state at J/W = 1 for the original necklace model. 

For general values of J/W, we solved (29) numerically. Fig. 1 shows the 
coupling dependence of T/v, for several values of N. On the horizontal axis, we 
have chosen the variable J / ( J + W), so as to be able to include the asymptotic 
regime where (30) becomes valid. (Since in the MFA, the coordination number 
z enters only through zW, we have simply taken z = 2). 

At J = 0, the system separates into N independent channels, disconnected 
from the S-spins. Each set of r-spins orders at the same mean-field tempera- 
ture W/2N. This point is clearly pathological; our interest lies in the behaviour 
at intermediate (and in any case, non-zero) values of J/W. 

First of all, one should note that, with the exception of N = 2, all curves 
show an overall rising tendency with increasing J/W: the asymptotic value 
(30) is lying significantly higher than the low- J value of T/v. Remarkably, the 
Kondo effect does not seem to compete with pseudopsin ordering but rather to 
assist it ! We can understand this tendency by remembering that the ordering 
is that of composite spins [most clearly seen in the large-J expression (32)] 
which arise from glueing the r-spins together. In the starting hamiltonian (3), 
there is no direct interchannel coupling: the effective ferromagnetic coupling 
between the channels is done by the Kondo coupling. So, in the multi-channel 
case, the Kondo coupling is effectively creating the objects which subsequently 
order. 

A subtler feature of our T N curves is a maximum at some J/W value 
which appears to increase with N. For a better view, we give a blow-up of 
the relevant region of the phase diagram in Fig. 2. The broad, relatively low 



21 



maximum is most discernable for N = 2, but it apparently exists (though 
quickly becoming quite inconspicuous) for any finite N. 

We have also done the numerical solution of the self-consistency equation 
at arbitary temperatures, i.e., at finite effective field strengths. We found that 
the temperature dependence of the order parameter is of the form usually 
found in similar mean-field solutions, so we renounce including them here. 

In our pseudospin model, the ordering which sets in at T N , is (depending on 
the sign of W), either a ferromagnetic, or a two-sublattice antiferromagnetic, 
ordering of the r-spins. Via the Kondo-coupling, this induces a similar, but 
oppositely polarized, ordering of the S'-spins. It depends on the nature of 
the underlying overscreened Kondo model, what the physical meaning of this 
order is. In Cox's ]17| model of U-compounds, it is primarily the ordering of 



electric quadrupole moments. The possibility that actual spin magnetization 
is a secondary order parameter arising from the mixing-in of higher-lying ionic 
levels, is an attractive possibility to account for the smallness of the ordered 
magnetic moment of some heavy fermion magnets [^4]]. It remains an open 
question whether quadrupolar Kondo effect can assist superconductivity . 



4 DISCUSSION AND SUMMARY 

The crucial question of the theory of Kondo lattices is to what extent can the 
results for a single Kondo impurity be taken as a guide for the behaviour of pe- 
riodic systems ? The answer must certainly depend on the size of the localized 
spin 5*, and the number of conduction electron channels N. We are uncer- 
tain about the answer even in the relatively simple cases of the underscreened 
(S > N/2), and exactly screened (S = N/2), Kondo lattices. The situation 
is quite obscure for the overscreened (S < N/2) Kondo lattice since it is dif- 
ficult to infer how the subtle, and puzzling, features revealed by the solution 
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of the single-ion multi-channel Kondo problem should manifest themselves in 
the lattice case. The available arguments suggest that switching on 

intersite interactions must have a drastic effect. 

Inspired by the fact that valuable insight into the behaviour of the S = 1/2, 
N = 1 Kondo lattice has been gained from studying the Kondo necklace 
model introduced by Doniach f25 |, we set out to define and investigate similar 
pseudospin models corresponding to a variety of Kondo lattices. Following our 
earlier study of underscreened Kondo necklaces [Q, here we introduced the N- 
channel, S = 1/2, overscreened necklace model (3). In the spirit of necklace 
models, the N screening channels of the underlying Kondo lattice problem 
are represented by N sets of pseudospins Tj. The low-lying electron-hole 
excitations of the Fermi sea, and hence the possibility of a "true" Kondo effect, 
are thereby lost. Of the conduction electrons, only their spin degrees of freedom 
are retained: the pseudospin model still incorporates the competing tendencies 
of a Kondo-like spin compensation, and magnetic ordering. Intersite coupling 
is mediated by the propagating character of the pseudospins, described by an 
x — y coupling term in (3). 

We studied the phase diagram of the model in the same kind of mean-field 
approximation (MFA) as that used by Doniach |23|. Setting up the self- 
consistency equation requires the diagonalization of the single-site effective 
hamiltonian (4) in a 2 JV+1 -dimensional Hilbert space: this is made possible 
by exploiting the invariance under permutations of the channel indices. We 
found that the low-temperature phase is always ordered; for W > 0, and 
bipartite lattices, this is just the Neel order of the r-spins, accompanied by 
the oppositely polarized ordering pattern of the £-spins. All the channels 
are equally polarized (5), corresponding to the well-known tendency that the 
Kondo-term mediates an effective ferromagnetic coupling between the different 
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channels. 

The nearly parallel alignment of the pseudospins results in a remarkable 
effect: for N > 3, the ordering tendency is actually enhanced with increasing 
Kondo coupling J ! This can be understood as arising from the fact that 
overscreening builds up composite pseudospin objects which are larger than 
the "naked" original spin S. Their coupling, mediated by the pseudospin x — y 
term, turns out to be sufficiently effective to increase the Neel temperature 
T N . 

In the terms of the model (3), the ordering is of the easy-plane type. This is 
most clearly born out by the exact form (32) of the effective hamiltonian which 
becomes valid in the limit of infinite Kondo coupling. Introducing the x — y 
form of the intersite r-coupling has destroyed the spin-rotational invariance 
of the underlying electronic hamiltonian, which is still shown by the isotropic 
Kondo-coupling. In this sense, the necklace models are in a universality class 
different from that of the Kondo lattices |2"7|] , and this has to be kept in mind 
when trying to transfer results from one class of models to the other. 

The low-temperature state of the model can be visualized as the order- 
ing of composite \r]\ = (N — l)/2 spins. This apparently corresponds to the 
strong-coupling behaviour of the overcompensated Kondo-site, and not to 
its intermediate-coupling fixed point with the strange, quantum-mechanical 
non-integer zero point degeneracy ]IUL [TTJ. The conjecture JE^I that the over- 
screened lattice might show a corresponding exotic phase is certainly exciting 
but this is not born out by our study of a drastically simplified model. On 



the other hand, we can relate our findings to independent arguments [pi] , |32 



suggesting that intersite interactions are likely to drive the system towards the 
strong-coupling regime. 

Noting that the multi-channel Kondo problem can be used to model a 
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variety of systems |T3|, [IT], [19], p0[ , the physical nature of the predicted order 
can also be different from case to case. Thinking of potential applications 
we have, however, primarily the suggested quadrupolar ordering of U-based 
systems in mind |17|, |24| . 

To summarize, we introduced the iV-channel Kondo necklace model for 
a preliminary study of the nature of the collective spin state in overscreened 
Kondo lattices. Our findings indicate that intersite interactions drastically 
change the state of the Kondo centres, lifting the ground state degeneracy, 
and inducing (pseudo)spin ordering at low temperatures. 
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Figure Captions 



Fig. 1 The dependence of the Neel temperature T N on the dimensionless 
Kondo coupling J/W, for iV=2, 3,4,5, and 10, and z — 2. Choosing the vari- 
able J / ( J + W) on the horizontal axis compresses the half-axis < J/W < oo 
into a finite interval. 

Fig. 2 Enlargement of part of Fig. 1 to show more clearly the maxima on 
the T/v vs J j (J + W) curves. The positions of the maxima are indicated by 
diamonds. The dotted line is a guide to the eye; it indicates that the maximum 
shifts to increasingly large values of J/W with increasing N . 
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